The impact of damaging epilepsy and cardiac genetic variant burden in sudden death in the young.

BACKGROUND: Sudden unexpected death in children is a tragic event. Understanding the genetics of sudden death in the young (SDY) enables family counseling and cascade screening. The objective of this study was to characterize genetic variation in an SDY cohort using whole genome sequencing. METHODS: The SDY Case Registry is a National Institutes of Health/Centers for Disease Control and Prevention surveillance effort to discern the prevalence, causes, and risk factors for SDY. The SDY Case Registry prospectively collected clinical data and DNA biospecimens from SDY cases < 20 years of age. SDY cases were collected from medical examiner and coroner offices spanning 13 US jurisdictions from 2015 to 2019. The cohort included 211 children (median age 0.33 year; range 0-20 years), determined to have died suddenly and unexpectedly and from whom DNA biospecimens for DNA extractions and next-of-kin consent were ascertained. A control cohort consisted of 211 randomly sampled, sex- and ancestry-matched individuals from the 1000 Genomes Project. Genetic variation was evaluated in epilepsy, cardiomyopathy, and arrhythmia genes in the SDY and control cohorts. American College of Medical Genetics/Genomics guidelines were used to classify variants as pathogenic or likely pathogenic. Additionally, pathogenic and likely pathogenic genetic variation was identified using a Bayesian-based artificial intelligence (AI) tool. RESULTS: The SDY cohort was 43% European, 29% African, 3% Asian, 16% Hispanic, and 9% with mixed ancestries and 39% female. Six percent of the cohort was found to harbor a pathogenic or likely pathogenic genetic variant in an epilepsy, cardiomyopathy, or arrhythmia gene. The genomes of SDY cases, but not controls, were enriched for rare, potentially damaging variants in epilepsy, cardiomyopathy, and arrhythmia-related genes. A greater number of rare epilepsy genetic variants correlated with younger age at death. CONCLUSIONS: While damaging cardiomyopathy and arrhythmia genes are recognized contributors to SDY, we also observed an enrichment in epilepsy-related genes in the SDY cohort and a correlation between rare epilepsy variation and younger age at death. These findings emphasize the importance of considering epilepsy genes when evaluating SDY.

Selection, optimization, and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse populations.

Polygenic risk scores (PRS) have improved in predictive performance supporting their use in clinical practice. Reduced predictive performance of PRS in diverse populations can exacerbate existing health disparities. The NHGRI-funded eMERGE Network is returning a PRS-based genome-informed risk assessment to 25,000 diverse adults and children. We assessed PRS performance, medical actionability, and potential clinical utility for 23 conditions. Standardized metrics were considered in the selection process with additional consideration given to strength of evidence in African and Hispanic populations. Ten conditions were selected with a range of high-risk thresholds: atrial fibrillation, breast cancer, chronic kidney disease, coronary heart disease, hypercholesterolemia, prostate cancer, asthma, type 1 diabetes, obesity, and type 2 diabetes. We developed a pipeline for clinical PRS implementation, used genetic ancestry to calibrate PRS mean and variance, created a framework for regulatory compliance, and developed a PRS clinical report. eMERGE's experience informs the infrastructure needed to implement PRS-based implementation in diverse clinical settings.

The impact of damaging epilepsy and cardiac genetic variant burden in sudden death in the young.

Background: Sudden unexpected death in children is a tragic event. Understanding the genetics of sudden death in the young (SDY) enables family counseling and cascade screening. The objective of this study was to characterize genetic variation in an SDY cohort using whole genome sequencing. Methods: The SDY Case Registry is a National Institutes of Health/Centers for Disease Control surveillance effort to discern the prevalence, causes, and risk factors for SDY. The SDY Case Registry prospectively collected clinical data and DNA biospecimens from SDY cases <20 years of age. SDY cases were collected from medical examiner and coroner offices spanning 13 US jurisdictions from 2015-2019. The cohort included 211 children (mean age 1 year; range 0-20 years), determined to have died suddenly and unexpectedly and in whom DNA biospecimens and next-of-kin consent were ascertained. A control cohort consisted of 211 randomly sampled, sex-and ancestry-matched individuals from the 1000 Genomes Project. Genetic variation was evaluated in epilepsy, cardiomyopathy and arrhythmia genes in the SDY and control cohorts. American College of Medical Genetics/Genomics guidelines were used to classify variants as pathogenic or likely pathogenic. Additionally, genetic variation predicted to be damaging was identified using a Bayesian-based artificial intelligence (AI) tool. Results: The SDY cohort was 42% European, 30% African, 17% Hispanic, and 11% with mixed ancestries, and 39% female. Six percent of the cohort was found to harbor a pathogenic or likely pathogenic genetic variant in an epilepsy, cardiomyopathy or arrhythmia gene. The genomes of SDY cases, but not controls, were enriched for rare, damaging variants in epilepsy, cardiomyopathy and arrhythmia-related genes. A greater number of rare epilepsy genetic variants correlated with younger age at death. Conclusions: While damaging cardiomyopathy and arrhythmia genes are recognized contributors to SDY, we also observed an enrichment in epilepsy-related genes in the SDY cohort, and a correlation between rare epilepsy variation and younger age at death. These findings emphasize the importance of considering epilepsy genes when evaluating SDY.

Opposing effects of genetic variation in MTCH2 for obesity versus heart failure.

Genetic variation in genes regulating metabolism may be advantageous in some settings but not others. The non-failing adult heart relies heavily on fatty acids as a fuel substrate and source of ATP. In contrast, the failing heart favors glucose as a fuel source. A bootstrap analysis for genes with deviant allele frequencies in cardiomyopathy cases versus controls identified the MTCH2 gene as having unusual variation. MTCH2 encodes an outer mitochondrial membrane protein, and prior genome-wide studies associated MTCH2 variants with body mass index, consistent with its role in metabolism. We identified the referent allele of rs1064608 (p.Pro290) as being overrepresented in cardiomyopathy cases compared to controls, and linkage disequilibrium analysis associated this variant with the MTCH2 cis eQTL rs10838738 and lower MTCH2 expression. To evaluate MTCH2, we knocked down Mtch in Drosophila heart tubes which produced a dilated and poorly functioning heart tube, reduced adiposity and shortened life span. Cardiac Mtch mutants generated more lactate at baseline, and they displayed impaired oxygen consumption in the presence of glucose but not palmitate. Treatment of cardiac Mtch mutants with dichloroacetate, a pyruvate dehydrogenase kinase inhibitor, reduced lactate and rescued lifespan. Deletion of MTCH2 in human cells similarly impaired oxygen consumption in the presence of glucose but not fatty acids. These data support a model in which MTCH2 reduction may be favorable when fatty acids are the major fuel source, favoring lean body mass. However, in settings like heart failure, where the heart shifts toward using more glucose, reduction of MTCH2 is maladaptive.

Bacterial Superinfection Pneumonia in Patients Mechanically Ventilated for COVID-19 Pneumonia.

Rationale: Current guidelines recommend patients with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pneumonia receive empirical antibiotics for suspected bacterial superinfection on the basis of weak evidence. Rates of ventilator-associated pneumonia (VAP) in clinical trials of patients with SARS-CoV-2 pneumonia are unexpectedly low. Objectives: We conducted an observational single-center study to determine the prevalence and etiology of bacterial superinfection at the time of initial intubation and the incidence and etiology of subsequent bacterial VAP in patients with severe SARS-CoV-2 pneumonia. Methods: Bronchoscopic BAL fluid samples from all patients with SARS-CoV-2 pneumonia requiring mechanical ventilation were analyzed using quantitative cultures and a multiplex PCR panel. Actual antibiotic use was compared with guideline-recommended therapy. Measurements and Main Results: We analyzed 386 BAL samples from 179 patients with SARS-CoV-2 pneumonia requiring mechanical ventilation. Bacterial superinfection within 48 hours of intubation was detected in 21% of patients. Seventy-two patients (44.4%) developed at least one VAP episode (VAP incidence rate = 45.2/1,000 ventilator days); 15 (20.8%) initial VAPs were caused by difficult-to-treat pathogens. The clinical criteria did not distinguish between patients with or without bacterial superinfection. BAL-based management was associated with significantly reduced antibiotic use compared with guideline recommendations. Conclusions: In patients with SARS-CoV-2 pneumonia requiring mechanical ventilation, bacterial superinfection at the time of intubation occurs in <25% of patients. Guideline-based empirical antibiotic management at the time of intubation results in antibiotic overuse. Bacterial VAP developed in 44% of patients and could not be accurately identified in the absence of microbiologic analysis of BAL fluid.

Genomic Autopsy of Sudden Deaths in Young Individuals.

Importance: Postmortem genetic testing of young individuals with sudden death has previously identified pathogenic gene variants. However, prior studies primarily considered highly penetrant monogenic variants, often without detailed decedent and family clinical information. Objective: To assess genotype and phenotype risk in a diverse cohort of young decedents with sudden death and their families. Design, Setting, and Participants: Pathological and whole-genome sequence analysis was conducted in a cohort referred from a national network of medical examiners. Cases were accrued prospectively from May 2015 to March 2019 across 24 US states. Analysis began September 2016 and ended November 2020. Exposures: Evaluation of autopsy and clinical data integrated with whole-genome sequence data and family member evaluation. Results: A total of 103 decedents (mean [SD] age at death, 23.7 [11.9] years; age range, 1-44 years), their surviving family members, and 140 sex- and genetic ancestry-matched controls were analyzed. Among 103 decedents, autopsy and clinical data review categorized 36 decedents with postmortem diagnoses, 23 decedents with findings of uncertain significance, and 44 with sudden unexplained death. Pathogenic/likely pathogenic (P/LP) genetic variants in arrhythmia or cardiomyopathy genes were identified in 13 decedents (12.6%). A multivariable analysis including decedent phenotype, ancestry, and sex demonstrated that younger decedents had a higher burden of P/LP variants and select variants of uncertain significance (effect size, -1.64; P = .001). These select, curated variants of uncertain significance in cardiac genes were more common in decedents than controls (83 of 103 decedents [86%] vs 100 of 140 controls [71%]; P = .005), and decedents harbored more rare cardiac variants than controls (2.3 variants per individual vs 1.8 in controls; P = .006). Genetic testing of 31 parent-decedent trios and 14 parent-decedent dyads revealed 8 transmitted P/LP variants and 1 de novo P/LP variant. Incomplete penetrance was present in 6 of 8 parents who transmitted a P/LP variant. Conclusions and Relevance: Whole-genome sequencing effectively identified P/LP variants in cases of sudden death in young individuals, implicating both arrhythmia and cardiomyopathy genes. Genomic analyses and familial phenotype association suggest potentially additive, oligogenic risk mechanisms for sudden death in this cohort.

Genomic Context Differs Between Human Dilated Cardiomyopathy and Hypertrophic Cardiomyopathy.

Background Inherited cardiomyopathies display variable penetrance and expression, and a component of phenotypic variation is genetically determined. To evaluate the genetic contribution to this variable expression, we compared protein coding variation in the genomes of those with hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM). Methods and Results Nonsynonymous single-nucleotide variants (nsSNVs) were ascertained using whole genome sequencing from familial cases of HCM (n=56) or DCM (n=70) and correlated with echocardiographic information. Focusing on nsSNVs in 102 genes linked to inherited cardiomyopathies, we correlated the number of nsSNVs per person with left ventricular measurements. Principal component analysis and generalized linear models were applied to identify the probability of cardiomyopathy type as it related to the number of nsSNVs in cardiomyopathy genes. The probability of having DCM significantly increased as the number of cardiomyopathy gene nsSNVs per person increased. The increase in nsSNVs in cardiomyopathy genes significantly associated with reduced left ventricular ejection fraction and increased left ventricular diameter for individuals carrying a DCM diagnosis, but not for those with HCM. Resampling was used to identify genes with aberrant cumulative allele frequencies, identifying potential modifier genes for cardiomyopathy. Conclusions Participants with DCM had more nsSNVs per person in cardiomyopathy genes than participants with HCM. The nsSNV burden in cardiomyopathy genes did not correlate with the probability or manifestation of left ventricular measures in HCM. These findings support the concept that increased variation in cardiomyopathy genes creates a genetic background that predisposes to DCM and increased disease severity.

Comparative Effectiveness of Surgical Approaches for Lung Cancer.

BACKGROUND: The magnitude of association and quality of evidence comparing surgical approaches for lung cancer resection has not been analyzed. This has resulted in conflicting information regarding the relative superiority of the different approaches and disparate opinions on the optimal surgical treatment. We reviewed and systematically analyzed all published data comparing near- (30-d) and long-term mortality for minimally invasive to open surgical approaches for lung cancer. METHODS: Comprehensive search of EMBASE, MEDLINE, and the Cochrane Library, from January 2009 to August 2019, was performed to identify the studies and those that passed bias assessment were included in the analysis utilizing propensity score matching techniques. Meta-analysis was performed using random-effects and fixed-effects models. Risk of bias was assessed via the Newcastle-Ottawa Scale and the ROBINS-I tool. The study was registered in PROSPERO (CRD42020150923) prior to analysis. RESULTS: Overall, 1382 publications were identified but 19 studies were included encompassing 47,054 patients after matching. Minimally invasive techniques were found to be superior with respect to near-term mortality in early and advanced-stage lung cancer (risk ratio 0.45, 95% confidence interval [CI] 0.21-0.95, I2 = 0%) as well as for elderly patients (odds ratio 0.45, 95% CI 0.31-0.65, I2 = 30%), but did not demonstrate benefit for high-risk patients (odds ratio 0.74, 95% CI 0.06-8.73, I2 = 78%). However, no difference was found in long-term survival. CONCLUSIONS: We performed the first systematic review and meta-analysis to compare surgical approaches for lung cancer which indicated that minimally invasive techniques may be superior to thoracotomy in near-term mortality, but there is no difference in long-term outcomes.

Pathogenic and Uncertain Genetic Variants Have Clinical Cardiac Correlates in Diverse Biobank Participants.

Background Genome sequencing coupled with electronic heath record data can uncover medically important genetic variation. Interpretation of rare genetic variation and its role in mediating cardiovascular phenotypes is confounded by variants of uncertain significance. Methods and Results We analyzed the whole genome sequence of 900 racially and ethnically diverse biobank participants selected from a single US center. Participants were equally divided among European, African, Hispanic, and mixed races/ethnicities. We evaluated the American College of Medical Genetics and Genomics medically actionable list of 59 genes, focusing on the cardiac genes. Variation was interpreted using the most recent reports in ClinVar, a database of medically relevant human variation. We identified 19 individuals with pathogenic or likely pathogenic variants in cardiac actionable genes (2%) and found evidence of related clinical correlates in the electronic health record. Participants of African ancestry, compared with those of European ancestry, had more variants of uncertain significance in the medically actionable genes including the 30 cardiac actionable genes, even when normalized to total variant count per person. Longitudinal measures of left ventricle size from ≈400 biobank participants (1723 patient-years) were correlated with genetic findings. The presence of ≥1 uncertain variant in the actionable cardiac genes and a cardiomyopathy diagnosis correlated with increased left ventricular internal diameter in diastole and in systole. In particular, MYBPC3 was identified as a gene with excess variants of uncertain significance. Conclusions These data indicate that a subset of uncertain genetic variants may confer risk and should not be considered benign.

Association of Cardiomyopathy With MYBPC3 D389V and MYBPC3Δ25bpIntronic Deletion in South Asian Descendants.

Importance: The genetic variant MYBPC3Δ25bp occurs in 4% of South Asian descendants, with an estimated 100 million carriers worldwide. MYBPC3 Δ25bp has been linked to cardiomyopathy and heart failure. However, the high prevalence of MYBPC3Δ25bp suggests that other stressors act in concert with MYBPC3Δ25bp. Objective: To determine whether there are additional genetic factors that contribute to the cardiomyopathic expression of MYBPC3Δ25bp. Design, Setting, andParticipants: South Asian individuals living in the United States were screened for MYBPC3Δ25bp, and a subgroup was clinically evaluated using electrocardiograms and echocardiograms at Loyola University, Chicago, Illinois, between January 2015 and July 2016. Main Outcomes and Measures: Next-generation sequencing of 174 cardiovascular disease genes was applied to identify additional modifying gene mutations and correlate genotype-phenotype parameters. Cardiomyocytes derived from human-induced pluripotent stem cells were established and examined to assess the role of MYBPC3Δ25bp. Results: In this genotype-phenotype study, individuals of South Asian descent living in the United States from both sexes (36.23% female) with a mean population age of 48.92 years (range, 18-84 years) were recruited. Genetic screening of 2401 US South Asian individuals found an MYBPC3Δ25bpcarrier frequency of 6%. A higher frequency of missense TTN variation was found in MYBPC3Δ25bp carriers compared with noncarriers, identifying distinct genetic backgrounds within the MYBPC3Δ25bp carrier group. Strikingly, 9.6% of MYBPC3Δ25bp carriers also had a novel MYBPC3 variant, D389V. Family studies documented D389V was in tandem on the same allele as MYBPC3Δ25bp, and D389V was only seen in the presence of MYBPC3Δ25bp. In contrast to MYBPC3Δ25bp, MYBPC3Δ25bp/D389V was associated with hyperdynamic left ventricular performance (mean [SEM] left ventricular ejection fraction, 66.7 [0.7%]; left ventricular fractional shortening, 36.6 [0.6%]; P < .03) and stem cell-derived cardiomyocytes exhibited cellular hypertrophy with abnormal Ca2+ transients. Conclusions and Relevance: MYBPC3Δ25bp/D389V is associated with hyperdynamic features, which are an early finding in hypertrophic cardiomyopathy and thought to reflect an unfavorable energetic state. These findings support that a subset of MYBPC3Δ25bp carriers, those with D389V, account for the increased risk attributed to MYBPC3Δ25bp.

Experimental Modeling Supports a Role for MyBP-HL as a Novel Myofilament Component in Arrhythmia and Dilated Cardiomyopathy.

BACKGROUND: Cardiomyopathy and arrhythmias are under significant genetic influence. Here, we studied a family with dilated cardiomyopathy and associated conduction system disease in whom prior clinical cardiac gene panel testing was unrevealing. METHODS: Whole-genome sequencing and induced pluripotent stem cells were used to examine a family with dilated cardiomyopathy and atrial and ventricular arrhythmias. We also characterized a mouse model with heterozygous and homozygous deletion of Mybphl. RESULTS: Whole-genome sequencing identified a premature stop codon, R255X, in the MYBPHL gene encoding MyBP-HL (myosin-binding protein-H like), a novel member of the myosin-binding protein family. MYBPHL was found to have high atrial expression with low ventricular expression. We determined that MyBP-HL protein was myofilament associated in the atria, and truncated MyBP-HL protein failed to incorporate into the myofilament. Human cell modeling demonstrated reduced expression from the mutant MYBPHL allele. Echocardiography of Mybphl heterozygous and null mouse hearts exhibited a 36% reduction in fractional shortening and an increased diastolic ventricular chamber size. Atria weight normalized to total heart weight was significantly increased in Mybphl heterozygous and null mice. Using a reporter system, we detected robust expression of Mybphl in the atria, and in discrete puncta throughout the right ventricular wall and septum, as well. Telemetric electrocardiogram recordings in Mybphl mice revealed cardiac conduction system abnormalities with aberrant atrioventricular conduction and an increased rate of arrhythmia in heterozygous and null mice. CONCLUSIONS: The findings of reduced ventricular function and conduction system defects in Mybphl mice support that MYBPHL truncations may increase risk for human arrhythmias and cardiomyopathy.

Multiscale Aspects of Generation of High-Gamma Activity during Seizures in Human Neocortex.

High-gamma (HG; 80-150 Hz) activity in macroscopic clinical records is considered a marker for critical brain regions involved in seizure initiation; it is correlated with pathological multiunit firing during neocortical seizures in the seizure core, an area identified by correlated multiunit spiking and low frequency seizure activity. However, the effects of the spatiotemporal dynamics of seizure on HG power generation are not well understood. Here, we studied HG generation and propagation, using a three-step, multiscale signal analysis and modeling approach. First, we analyzed concurrent neuronal and microscopic network HG activity in neocortical slices from seven intractable epilepsy patients. We found HG activity in these networks, especially when neurons displayed paroxysmal depolarization shifts and network activity was highly synchronized. Second, we examined HG activity acquired with microelectrode arrays recorded during human seizures (n = 8). We confirmed the presence of synchronized HG power across microelectrode records and the macroscale, both specifically associated with the core region of the seizure. Third, we used volume conduction-based modeling to relate HG activity and network synchrony at different network scales. We showed that local HG oscillations require high levels of synchrony to cross scales, and that this requirement is met at the microscopic scale, but not within macroscopic networks. Instead, we present evidence that HG power at the macroscale may result from harmonics of ongoing seizure activity. Ictal HG power marks the seizure core, but the generating mechanism can differ across spatial scales.

Network burst activity in hippocampal neuronal cultures: the role of synaptic and intrinsic currents.

The goal of this work was to define the contributions of intrinsic and synaptic mechanisms toward spontaneous network-wide bursting activity, observed in dissociated rat hippocampal cell cultures. This network behavior is typically characterized by short-duration bursts, separated by order of magnitude longer interburst intervals. We hypothesize that while short-timescale synaptic processes modulate spectro-temporal intraburst properties and network-wide burst propagation, much longer timescales of intrinsic membrane properties such as persistent sodium (Nap) currents govern burst onset during interburst intervals. To test this, we used synaptic receptor antagonists picrotoxin, 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX), and 3-(2-carboxypiperazine-4-yl)propyl-1-phosphonate (CPP) to selectively block GABAA, AMPA, and NMDA receptors and riluzole to selectively block Nap channels. We systematically compared intracellular activity (recorded with patch clamp) and network activity (recorded with multielectrode arrays) in eight different synaptic connectivity conditions: GABAA + NMDA + AMPA, NMDA + AMPA, GABAA + AMPA, GABAA + NMDA, AMPA, NMDA, GABAA, and all receptors blocked. Furthermore, we used mixed-effects modeling to quantify the aforementioned independent and interactive synaptic receptor contributions toward spectro-temporal burst properties including intraburst spike rate, burst activity index, burst duration, power in the local field potential, network connectivity, and transmission delays. We found that blocking intrinsic Nap currents completely abolished bursting activity, demonstrating their critical role in burst onset within the network. On the other hand, blocking different combinations of synaptic receptors revealed that spectro-temporal burst properties are uniquely associated with synaptic functionality and that excitatory connectivity is necessary for the presence of network-wide bursting. In addition to confirming the critical contribution of direct excitatory effects, mixed-effects modeling also revealed distinct combined (nonlinear) contributions of excitatory and inhibitory synaptic activity to network bursting properties.

Zodiac: A Comprehensive Depiction of Genetic Interactions in Cancer by Integrating TCGA Data.

BACKGROUND: Genetic interactions play a critical role in cancer development. Existing knowledge about cancer genetic interactions is incomplete, especially lacking evidences derived from large-scale cancer genomics data. The Cancer Genome Atlas (TCGA) produces multimodal measurements across genomics and features of thousands of tumors, which provide an unprecedented opportunity to investigate the interplays of genes in cancer. METHODS: We introduce Zodiac, a computational tool and resource to integrate existing knowledge about cancer genetic interactions with new information contained in TCGA data. It is an evolution of existing knowledge by treating it as a prior graph, integrating it with a likelihood model derived by Bayesian graphical model based on TCGA data, and producing a posterior graph as updated and data-enhanced knowledge. In short, Zodiac realizes "Prior interaction map + TCGA data → Posterior interaction map." RESULTS: Zodiac provides molecular interactions for about 200 million pairs of genes. All the results are generated from a big-data analysis and organized into a comprehensive database allowing customized search. In addition, Zodiac provides data processing and analysis tools that allow users to customize the prior networks and update the genetic pathways of their interest. Zodiac is publicly available at www.compgenome.org/ZODIAC. CONCLUSIONS: Zodiac recapitulates and extends existing knowledge of molecular interactions in cancer. It can be used to explore novel gene-gene interactions, transcriptional regulation, and other types of molecular interplays in cancer.

Estimating screening-mammography receiver operating characteristic (ROC) curves from stratified random samples of screening mammograms: a simulation study.

RATIONALE AND OBJECTIVES: To evaluate stratified random sampling (SRS) of screening mammograms by (1) Breast Imaging Reporting and Data System (BI-RADS) assessment categories, and (2) the presence of breast cancer in mammograms, for estimation of screening-mammography receiver operating characteristic (ROC) curves in retrospective observer studies. MATERIALS AND METHODS: We compared observer study case sets constructed by (1) random sampling (RS); (2) SRS with proportional allocation (SRS-P) with BI-RADS 1 and 2 noncancer cases accounting for 90.6% of all noncancer cases; (3) SRS with disproportional allocation (SRS-D) with BI-RADS 1 and 2 noncancer cases accounting for 10%-80%; and (4) SRS-D and multiple imputation (SRS-D + MI) with missing BI-RADS 1 and 2 noncancer cases imputed to recover the 90.6% proportion. Monte Carlo simulated case sets were drawn from a large case population modeled after published Digital Mammography Imaging Screening Trial data. We compared the bias, root-mean-square error, and coverage of 95% confidence intervals of area under the ROC curve (AUC) estimates from the sampling methods (200-2000 cases, of which 25% were cancer cases) versus from the large case population. RESULTS: AUC estimates were unbiased from RS, SRS-P, and SRS-D + MI, but biased from SRS-D. AUC estimates from SRS-P and SRS-D + MI had 10% smaller root-mean-square error than RS. CONCLUSIONS: Both SRS-P and SRS-D + MI can be used to obtain unbiased and 10% more efficient estimate of screening-mammography ROC curves.

Targeted analysis of whole genome sequence data to diagnose genetic cardiomyopathy.

BACKGROUND: Cardiomyopathy is highly heritable but genetically diverse. At present, genetic testing for cardiomyopathy uses targeted sequencing to simultaneously assess the coding regions of >50 genes. New genes are routinely added to panels to improve the diagnostic yield. With the anticipated $1000 genome, it is expected that genetic testing will shift toward comprehensive genome sequencing accompanied by targeted gene analysis. Therefore, we assessed the reliability of whole genome sequencing and targeted analysis to identify cardiomyopathy variants in 11 subjects with cardiomyopathy. METHODS AND RESULTS: Whole genome sequencing with an average of 37× coverage was combined with targeted analysis focused on 204 genes linked to cardiomyopathy. Genetic variants were scored using multiple prediction algorithms combined with frequency data from public databases. This pipeline yielded 1 to 14 potentially pathogenic variants per individual. Variants were further analyzed using clinical criteria and segregation analysis, where available. Three of 3 previously identified primary mutations were detected by this analysis. In 6 subjects for whom the primary mutation was previously unknown, we identified mutations that segregated with disease, had clinical correlates, and had additional pathological correlation to provide evidence for causality. For 2 subjects with previously known primary mutations, we identified additional variants that may act as modifiers of disease severity. In total, we identified the likely pathological mutation in 9 of 11 (82%) subjects. CONCLUSIONS: These pilot data demonstrate that ≈30 to 40× coverage whole genome sequencing combined with targeted analysis is feasible and sensitive to identify rare variants in cardiomyopathy-associated genes.

Supercomputing for the parallelization of whole genome analysis.

MOTIVATION: The declining cost of generating DNA sequence is promoting an increase in whole genome sequencing, especially as applied to the human genome. Whole genome analysis requires the alignment and comparison of raw sequence data, and results in a computational bottleneck because of limited ability to analyze multiple genomes simultaneously. RESULTS: We now adapted a Cray XE6 supercomputer to achieve the parallelization required for concurrent multiple genome analysis. This approach not only markedly speeds computational time but also results in increased usable sequence per genome. Relying on publically available software, the Cray XE6 has the capacity to align and call variants on 240 whole genomes in ∼50 h. Multisample variant calling is also accelerated. AVAILABILITY AND IMPLEMENTATION: The MegaSeq workflow is designed to harness the size and memory of the Cray XE6, housed at Argonne National Laboratory, for whole genome analysis in a platform designed to better match current and emerging sequencing volume.

Estimating sensitivity and specificity for technology assessment based on observer studies.

RATIONALE AND OBJECTIVES: The goal of this study was to determine the accuracy and precision of using scores from a receiver operating characteristic rating scale to estimate sensitivity and specificity. MATERIALS AND METHODS: We used data collected in a previous study that measured the improvements in radiologists' ability to classify mammographic microcalcification clusters as benign or malignant with and without the use of a computer-aided diagnosis scheme. Sensitivity and specificity were estimated from the rating data from a question that directly asked the radiologists their biopsy recommendations, which was used as the "truth," because it is the actual recall decision, thus it is their subjective truth. By thresholding the rating data, sensitivity and specificity were estimated for different threshold values. RESULTS: Because of interreader and intrareader variability, estimated sensitivity and specificity values for individual readers could be as much as 100% in error when using rating data compared to using the biopsy recommendation data. When pooled together, the estimates using thresholding the rating data were in good agreement with sensitivity and specificity estimated from the recommendation data. However, the statistical power of the rating data estimates was lower. CONCLUSIONS: By simply asking the observer his or her explicit recommendation (eg, biopsy or no biopsy), sensitivity and specificity can be measured directly, giving a more accurate description of empirical variability and the power of the study can be maximized.

Interreader scoring variability in an observer study using dual-modality imaging for breast cancer detection in women with dense breasts.

RATIONALE AND OBJECTIVES: To evaluate variability in the clinical assessment of breast images, we evaluated scoring behavior of radiologists in a retrospective reader study combining x-ray mammography (XRM) and three-dimensional automated breast ultrasound (ABUS) for breast cancer detection in women with dense breasts. METHODS: The study involved 17 breast radiologists in a sequential study design with readers first interpreting XRM-alone followed by an interpretation of combined XRM + ABUS. Each interpretation included a forced Breast Imaging Reporting and Data System scale and a likelihood that the woman had breast cancer. The analysis included 164 asymptomatic patients, including 31 breast cancer patients, with dense breasts and a negative screening XRM. Of interest were interreader scoring variability for XRM-alone, XRM + ABUS, and the sequential effect. In addition, a simulated double reading by pairs of readers of XRM + ABUS was investigated. Performance analysis included receiver operating characteristic analysis, percentile analysis, and κ statistics. Bootstrapping was used to determine statistical significance. RESULTS: The median change in area under the receiver operating characteristic curve after ABUS interpretation was 0.12 (range 0.04-0.19). Reader agreement was fair with the median interreader κ being 0.26 (0.05-0.48) for XRM-alone and 0.34 (0.11-0.55) for XRM + ABUS (95% confidence interval for the difference in κ, 0.06-0.11). Simulated double reading of XRM + ABUS demonstrated tradeoffs in sensitivity and specificity, but conservative simulated double reading resulted in a significant improvement in both sensitivity (16.7%) and specificity (7.6%) with respect to XRM-alone. CONCLUSION: A modest, but statistically significant, increase in interreader agreement was observed after interpretation of ABUS.

Receiver-operating characteristic curves for somatic cell scores and California mastitis test in Valle del Belice dairy sheep.

Using receiver-operating characteristic (ROC) curve methodology this study was designed to assess the diagnostic effectiveness of somatic cell count (SCC) and the California mastitis test (CMT) in Valle del Belice sheep, and to propose and evaluate threshold values for those tests that would optimally discriminate between healthy and infected udders. Milk samples (n=1357) were collected from 684 sheep in four flocks. The prevalence of infection, as determined by positive bacterial culture was 0.36, 87.7% of which were minor and 12.3% major pathogens. Of the culture negative samples, 83.7% had an SCC<500,000/mL and 97.4% had <1,000,000cells/mL. When the associations between SC score (SCS) and whole sample status (culture negative vs. infected), minor pathogen status (culture negative vs. infected with minor pathogens), major pathogen status (culture negative vs. infected with major pathogens), and CMT results were evaluated, the estimated area under the ROC curve was greater for glands infected with major compared to minor pathogens (0.88 vs. 0.73), whereas the area under the curve considering all pathogens was similar to the one for minor pathogens (0.75). The estimated optimal thresholds were 3.00 (CMT), 2.81 (SCS for the whole sample), 2.81 (SCS for minor pathogens), and 3.33 (SCS for major pathogens). These correctly classified, respectively, 69.0%, 73.5%, 72.6% and 91.0% of infected udders in the samples. The CMT appeared only to discriminate udders infected with major pathogens. In this population, SCS appeared to be the best indirect test of the bacteriological status of the udder.